REGULARITY AND CONVERGENCE RATES FOR THE LYAPUNOV 
EXPONENTS OF LINEAR CO-CYCLES 
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Abstract. We consider co-cycles A : TxK d e (x, v) h> (x+oj,A(x,E)v) with 
<o Diophantine, K = R or K = C. We assume that A : T x £ -> GL(d,K) 
is continuous, depends analytically on x 6 T and is Holder in E e £, 
where £ is a compact metric space. It is shown that if all Lyapunov 
exponents are distinct at one point E € £, then they remain distinct 
near E. Moreover, they depend in a Holder fashion on E e B for any ball 
B c £ where they are distinct. Similar results, with a weaker modulus of 
continuity, hold for higher-dimensional tori TP' with a Diophantine shift. 
We also derive optimal statements about the rate of convergence of the 
finite-scale Lyapunov exponents to their infinite-scale counterparts. A 
key ingredient in our arguments is the Avalanche Principle, a deterministic 
statement about long finite products of invertible matrices, which goes 
back to work of Michael Goldstein and the author. We also discuss 
applications of our techniques to products of random matrices. 



1. Introduction 

Let (X, £, m; 6) be an ergodic space, where m is an invariant measure 
under the ergodic transformation 6. Let A : X — > GL(d, K) where K — R or 
K — C, be a measurable transformation such that 

log||A(x)||, logllA^WIIe L\X,m) 

Define the co-cycle as 

A : X x K d 6 (x, i7) h+ (6x,A(x)v) 
Powers of A lead to matrix products 

4" } :=flA(9'x). 



The classical Oseledec theorem |Ose|, also known as multiplicative ergodic 
theorem, guarantees a filtration by linear subspaces, i.e., for a.e. x e X there 
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exist linear subspaces 

K d = V r x DV r ~ 1 D...DV 1 x DV x = {0} (1) 

and a sequence of positive numbers \i r > > ... > pi associated with 
this filtration, which determine the Lyapunov exponents as follows. One 

has that dim V x is a.e. constant, and furthermore, the (set) complement of 

V x 1 in V x is characterized by the property that 

ilogHA^H^log^ veVlwt 1 
as n — > oo. The Lyapunov exponents 

Ai > A 2 > • • • > A d -l > A rf (2) 



are precisely the numbers {log y-M =r counted with multiplicity given by 

dim V x - dim V x 
||Rag]lRueHLedl. 



dim V x - dim V x 1 . Finally, we remark that A(x)v{ = V 3 dx for a.e. x. See also 



The Lyapunov exponents are of fundamental importance in dynamics 
and have been widely studied, see the encyclopedic treatment in [ Bar Pes |. 
For more recent research literature see for example HAviVialt |AviVia2[ 
IBonViat IBGMVL Results for co-cycles with shift base (but in a very dif- 
ferent spirit from this paper) can be found in [Avi, AviKri] FayKri|lKri||. 



2. Statement of the main results 

We mainly consider co-cycles whose base dynamics exhibits very weak 
mixing properties. Let A : T — > GL(d, K) be continuous, and define the 
co-cycle 

1x^6 (x, v) m> (x + a),A(x)v) (3) 
where co is Diophantine, which will mean that 

c(co) 

\\na)\\ > -r -rz V n > 2 (4) 

with a > 1 arbitrary but fixed. Lebesgue almost every co € (0, 1) satisfy this 
condition. The purpose of this note is to point out a mild regularity result 
that the Lyapunov exponents exhibit as functions of a parameter which A 
is assumed to depend on. To be precise, we prove the following theorem 
where (£, b) is a compact metric space. The study of Lyapunov exponents 
(especially, conditions which ensure positivity of the top exponent) of linear 
co-cycles with deterministic base dynamics goes back to the seminal work 
by M. Herman ffHer| . 

Theorem 1. Suppose A : T x £ — » GL(d, K) is continuous (K - R, C), and 
analytic as a function x i-> A(x, E) uniformly in E e S. Furthermore, suppose 
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E i — > A(x, E) is Holder continuous, uniformly in x eT. Assume that the Eyapunov 
exponents satisfy the gap condition 

Aj(E) - A j+1 (E) > k > V E e 6, V 1 < ;' < d (5) 

Then all Ay(E) are Holder continuous as a function ofE e S. Moreover, if (O holds 
at some point Eo e £, then each A,(E) is Holder continuous locally around Eo. In 
other words, if all exponents are distinct at Eo, then they are all Holder continuous 
locally around Eq, and therefore also remain distinct near Eq. 

The proof is based on an adaptation of the technique developed jointly 
with Michael Goldstein in [GolSchl] for SL(2, R) Schrodinger co-cycles. In 
fact, [GolSchl J dealt with Schrodinger co-cycles for which E plays the role of 
the "energy" (spectral parameter), and one needed to impose the condition 
of positive Lyapunov exponents. In Theorem [TJ this positivity is exactly 
mirrored by the gap condition. For Schrodinger co-cycles Holder regularity 
of the exponents is best possible; this follows from the connection between 
the Lyapunov exponent and the integrated density of states as expressed 
by the Thouless formula, see IBou"T | . 

The Holder C a class which our argument produces for the Lyapunov 
exponents depends on the size of the gaps, i.e., a depends on k in ((5]). Note 
that for SL(2, R) co-cycles for which the larger exponent is positive, one can 
refine the argument in [GolSchl J in such a way that a does not depend on k, 
see |Bou2[ Chapter 8]. The point here is that for such SL(2, R) co-cycles the 
exponents are A(E) > > -A(E) whence one can take k > A. This property 
allows one to obtain uniform control on the Holder exponent a even as 
k — > 0. For GE(d, K) with d > 2 there is no relation between the size of 
the gaps and the individual exponents, and it therefore remains an open 
question as to how a may deteriorate as k — > 0, let alone what happens 
when the gap condition is violated. For example, one might expect that 
the top exponent remains Holder continuous provided 10 holds with j = 1. 
However, the argument presented in this paper breaks down completely if, 
say, Ai > A2 = A3. Two more comments are in order: 

• In addition to the regularity of the exponents, we also obtain an 
estimate on the rate of convergence of the finite-scale Lyapunov 
exponents. To be specific, we show in Section [6] that for all n > 1 

|Ai»(E) - Ay(E)| < - V£e£ (6) 

where Aj itl are the finite volume exponents. Moreover, one has the 
following dichotomy: (i) either (O is optimal along a sequence of 
n — > 00 or (ii) the convergence in (O is exponentially fast. 

• It is of course natural to inquire about other types of dynamics 
with weak mixing properties, such as shifts x h-» x + co on higher- 
dimensional tori TP, as well as shifts on T for which co satisfies 
weaker Diophantine conditions than Q). In these cases one has 
a somewhat weaker result than Theorem [TJ While the Lyapunov 
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exponents are again continuous, they are not seen to be Holder by 
the technology used here. Rather, one obtains regularity of the form 

|A/(E) - Aj(E')\ < C exp ( - c | log b(E, E')f) (7) 

with some < a < 1 and c > 0. See Section[7l 
We remark that 10 implies for all 1 < j < d 

dim V 1 X (E) - dim V j ~ l {E) = 1 for a.e. x £ T, V E e S (8) 

where {V ! x (E)}j denotes the Oseledec filtration in dependence of E. 

The technique developed in [GolSchl ) is based on the following two main 
ingredients: 

• The Avalanche Principle (AP) for 2 x 2 matrices. This is a deterministic 
statement. 

• A Earge Deviation Theorem (EDT) for the matrices A%. These are 
quantitative versions of the Furstenberg-Kesten theorem, or of King- 
man's sub-additive ergodic theorem. The LDTs are analytical state- 
ments, and they depend on the structure of the co-cycles and the 
dynamics on the base. 

The LDT theorem from BGolSchlL see also (BouGol| and |Bou2|, applies 
easily to higher-dimensional co-cycles. On the other hand, the AP does not 
carry over from the 2x2 case without changes. In fact, the gap condition is 
intimately tied to the formulation of the AP as it appears in the following 
section. In the final section of this paper we apply the main strategy to the 
products of random matrices and derive the same type of conclusions as 
we did for the shift dynamics. 

This note is organized as follows. In Section|3]we discuss the deterministic 
AP. In Section 0] we discuss the large deviation estimates for shifts on T 
under the Diophantine condition (H)). Here the analytic dependence of A on 
x is important. The proof of Theorem[T]is presented in Section[5j and sharp 
convergence rates for the Lyapunov exponents are presented in Section HI 
The final two sections are very sketchy, and discuss other types of dynamics 
on the base. While Section[7] deals with more general shifts, Section[8]covers 
products or random matrices. 

3. The Avalanche Principle 

The point of Lemma Q] below, which is a completely deterministic result, 
is that it effectively allows one to linearize long (non-commuting) matrix 
products A n A n _i . . . A^A\ . One cannot expect such a result for general 
products, even if each factor Aj has large norm. The mechanism which 
underlies this linearization is as follows: 

• we assume that each matrix has a dominant simple singular value, 
see (HI). 

• no two adjacent pairs Aj + \Aj cancel, in the sense that the dominant 
stretching of Aj is annihilated by Aj+i, see (flu). 
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Appropriate quantitative formulations of these properties guarantee that 
the whole product A n ...A\ is "to leading order" nothing other than the 
composition of the 1-dimensional dominant actions, see ((TTjl . This is what 
we mean by linearization. 

Lemma 1. Let {Ay}" =1 c GL(d,K) satisfy the following properties: for each 1 < 
j <n there exists a 1-dimensional subspace Sj c K d such that 

\AjV\ = \\Aj\\ \v\ Vv 6 Sj 

\Ajio\ < a.j\w\ Vw e Sj- (9) 
\\Aj\\ > ctj\i, p > 16n 2 

In addition, we assume that 

\\A j+1 \\ \\Aj\\ < pkA j+ iAj\\ V 1 < / < n. (10) 

Then 

n-l n-1 

log||A n -...-Ai|| + y log||Ay||-Y log||A ;+ iA ; || <C^- (11) 

;=2 ;=1 ^' 

for some absolute constant C. 



Proof. We denote by rij the orthogonal projection onto Sj, and by nj- that 
onto S\ Define 7y :- AjSj and 

Ay := —j- — , veTj\{0\ (12) 

Since Aj is invertible, we see that Tj is a line. By a polar decomposition, 
each Aj = UjDjUj where Dy is diagonal, and !iy, ITy are unitary. Moreover, 

Dy = diag(ayi, Oj /2 , . . . ,o ]A ), o j/k > V 1 < k < d 



and 



Then 



whence 



maxoyj- = 1 1 Ay 1 1 = ay^., maxayjt < a 



\\A j+1 Aj\\ < Aj\\A j+1 \\ \\Aj\\ + aj\\Aj +1 \\ + a j+1 \\Aj\ 
\\A j+1 Aj\\ > Aj\\A j+1 \\ \\Aj\\ - a j+1 \\Aj\\ 



i i l|A;+iA;|| i 

Pj - 2,- < Ay < Pj + ^ py := > ^ (13) 



In particular, 

Ay > - 2p~ X 
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SetM„ := nL„A ; -- Then for any u e K d 



£!,...,£„ 



where each £ ; is either the empty symbol, or equals _L. Consider the contri- 
bution to ([14]) where all £/ equal the empty symbol: 



\A n n n A n -iTl n -l O . . . O 7I 3 A 2 7l2Ai7ZiW| 

)3-l )3 

=i7HMij] A /ri 1 ^ 11 

7=1 7=1 



(15) 



Now let 7 := {;' € [1, n] | £ ; - =_L}, 7 C := [1, n] \ J, and estimate 

\A n n E «A n - X n e *z\ 0...0 A^A 1 n^A 1 n\ l u\ 

<\u\]\\\Aj\\Y[a h \[ Aj (16) 

;'el c kel {j+l,j)eFxF 

Comparing ((16)) with ((15)) yields 

(n^n^rn^n^ n ^ 

;=1 ;'=1 ;'eF fceZ (;+l,/)e7 c x/ c 

<^|( n A^ 1 <^(^4_ 2/ u- i r 2lfl <4'V-" 

/eZu(I-l) 

provided > 4/j -1 which is the case for \i > 8, say. Next, we note that 

(^) 4 V^ = (1 + 4/-i"5)" - 1 < 1271/U-5 

provided 4n/u~3 < 1, or /u > 16n 2 . Here we used that 

(1 + x)" - 1 = exp(n log(l + x)) - 1 < exp(nx) - 1 < 3nx if nx < 1 
In conclusion, 

(17) 

n-X n y ' 



> 



;'=1 ;=1 f=l 
n-1 n 

(l-12n/i-J)J]A 7 -J]|lAy|| 
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as well as 



n-l 



\\M n \\ < (1 + llnprl) Yl A; [] ||A ; || (18) 
;'=1 y=i 



Next, we note that 



n-l n-l 



(19) 



7=1 7=1 
Combining (O with Q7J) and I® implies pi). □ 

Note that Lemma [TJ remains unchanged if we replace any Ay with a 
nonzero multiple of itself. Also, for the scalar case ((TT]) is identically zero. 
The fact that Sj are lines, rather than planes or higher-dimensional sub- 
spaces, played a crucial role in the proof. Indeed, for higher-dimensional 
spaces we would need to replace l|T2)l with 

T \Uj+iV\ 

Aj := sup — — — 

veTj\{0} \ v \ 

Then ( fl3l) remains valid with A; instead of Aj. However, in the estimation 
of the dominant branch of ((14)) these quantities are not sufficient. Instead, 
we could only use the lower bounds 

A.:= inf !^ 

— J veTj\{0] \v\ 

and it is not possible to conclude as we did for dim(*Sy) = 1. Note that for 
the case of SL(2, R) matrices as in BGolSchlH the line-condition is automatic. 
To be precise, we simply take aj - so that condition I© becomes 

\\Aj\\>4t[i. 

To illustrate this point further, consider orthogonal projections Pj in JR 3 
or rank 1. Let 

Aj = Pj + £ (1 - Pj) (20) 

where e > is small but fixed. This ensures that Aj € GL(3, R). Then 
\\Aj\\ = 1, and ||A/+iA,|| = | cos 0j\ as e — > where 6j is the angle between the 
range of Pj and that of Pj+\ . Moreover, 

n-l 

||A H -....Ai||-» JJlcosS;-! 

7=1 
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as £ — > whence the expression in ((TTjl reads 

n— 1 n-1 

log \\A n ■ . . . ■ Ai| | + !og - Yj 1o S 

7=2 ;=1 
n-1 

= log \\A n AiH - £ log -> as e -> 

pi 

Lemma [I] gives the precise conditions in terms of (flO|) and I© such that ((TTjl 
holds for non-zero a in ((20]). Note that the former condition simply reads 

| COS 6j\ > |U _ 5. 

On the other hand, define now 

Aj := 1 - P ; - + ePy 

In that case = 1 for any < e < 1 since 2-planes have a non-zero 

intersection in R 3 . It follows that (flOT > trivially holds, as does (|9]) at least 
when e — > 0. However, ([TT]) becomes 

||A B -...-Ai||>exp(-C-^=) 

for any t > which is absurd since the left-hand side can be zero in the 
limit £ — > 0, and we may also send y. — » oo as £ — > 0. 



4. The Large Deviation Estimates 

Next, we state the LDTs for the co-cycles, see j Gol Schl , Theorem 3.8, 
Lemma 4.1], and |Bou2, Chapter 5]. 

Lemma 2. Under the Diophantine condition on co, see (@]), one has for any 5 > 
|{x 6 T | In" 1 log } (E)|| - Ai /n (£)| > d) < exp(-c6 2 n + C(log nf) (21) 
where c, C, b are positive constants depending on A, co, and n>2. Here 

Ax, n (E):= [ n- l \oz\\Af\E)\\dx (22) 
Jt 

for alln>l,Ee 6. 

Proof. The argument proceeds in two steps, the first step being analytical. 
The map x h-» n -1 log||A^(E)|| extends to a neighborhood J?l of T in the 
complex plane as a subharmonic function, uniformly in E e S. Moreover, we 
have 

supn^logHA^FOIUM VEe£ (23) 



r 

Jt 
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for some M = M(A, £) as well as 
n^logllAf^EJUdx 

> f (£fn)- 1 log(c(rf)|det^ n) (£)|)^ (24) 
Jt 

- I d" 1 log | det A(x, E)\ dx + (dn)" 1 log c(d) > -M VEefi 
Jt 

Properties d23l and ((24]) ensure that Theorem 3.8 in [GolSchl] applies to the 
subharmonic function 

u n (z,E) -(Mn)" 1 \og\\Af\E)\\ 
uniformly in n. The conclusion is as follows: for all 5 > 0, and any m > 1, 



(25) 



|{x e T | \rri~ 1 ^ «„(x + /a;, E) - (u n (-, E))\ > d] 

< exp(-c6m + C(logm) fc ) 

where (/) := L /(x) is the mean, and c, C, only depend on A, co (and not 
on m,8,n,E). 

The second step is to note the almost invariance of the co-cycles: 
sup In" 1 log |K1(E)H - n- 1 log ||4 n) (E)|| 



xeT 
< 



n-VsuplogllAfoE)!! + suplog||A- 1 (x / E)||) = - 

xeT xeT n 



Upon iteration, we have 

Cfc 

sup |w n (x + fay, E) - u„(x, E)| < — (26) 
xeT n 

To obtain (|2T|) we now set m := 6n/C in d25l) and conclude via d26t . □ 
By general principles, 

Ai /f! (E) — > Ai(E) as n — > oo 

as a monotone decreasing sequence. Indeed, it follows from the co-cycle 
relation 



that 



An+m) _ An) Am) 
A* - A x+mco ° At 



log||Ar m) ll<log|l41j|+log||4 m) l 



Averaging in x implies (n + m)A\ in+m < nX\ iU + mA\ im , whence the claim. 

We now turn to the exterior powers of A { "\ which we denote by APA ( X '. 
On a general vector space V with a linear operator L : V — > V one defines 

(A p L)/(£i, . . := f{Vi x , Vl v ) V ^ e V* 
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where / is an alternating linear form on (V*) p . We have the following 
more general version of the LDT, which allows us to control the Lyapunov 
exponents Ay with j > 2. 

Lemma 3. Under a suitable Diophantine condition on co one has for any 5 > 

v 

\[x e T | log||AP4' !) (£)H - Yj A '> (E) I > 6 ) < e M~ c ^ n + C(logn) 6 ) 

i=\ 

where c, C, b are positive constants depending on A, co, and n > 2. The A i/n are 
defined inductively via 



r 1 p 

-log||AMf (E)||dx = V A iin (E) 



i=l 

with 1 < p < d. 

Proof. The proof is essentially the same as that of the previous lemma; 
indeed, we retain the analyticity of 



on the annulus J?l uniformly in E € S. This allows us to apply the subhar- 
monic machinery as before. Moreover, we have the product structure 



AMf(E) = Yl A p A(x + jco, E) 



]=n 



which ensures the almost invariance (|26j) (note that A? A remains invertible). 
So the argument leading to the previous lemma applies. □ 

Write (dropping E for simplicity) A { " ] = U { " ] o D { " ] o U ( " ] , where U { "\ U { " ] 

(n) 

are unitary, and D x is diagonal with entries 

O hn (x) > a 2 , n (x) > . . . > Od-l,n( X ) ^ °d,n( x ) > 

It follows that APA^ = A P G^ o AW^ o AW™, whence 
f - log || AM$° || dx = f - log 1 1 APD™ 1 1 dx 

P 1 f P 
= y*- I log Oi, n {x)dx = V A,, M 

since A^U is an isometry if U is. Consequently, 

r !/n (x) = A;,„ V 1 < l < d (27) 



^X ioga - ( 
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In addition, we note that as decreasing limits of continuous functions, the 
sums 

V 

£/;(E), l<p<d 
i=l 

are upper-semicontinuous. But more regularity is harder to obtain. 

5. The proof of Theorem [T] 

Proof of Theorem\T\ Normalizing, we may assume that A : T — » SL(d,K). 
This implies 

d 

£a 7 -(E) = (28) 
;=i 

For simplicity, we consider the case of the two largest Lyapunov exponents 
first, and our goal is to show that the largest Lyapunov exponent is Holder. 
In other words, / = 1 in ©. Fix Eq e S. We may choose uq such that 

|A //n (£ ) - Ai(£o)l < 6o V n > n , 1 < i < d (29) 

where §o > is a small constant which will be specified later. By (|27| 

- J log||A^(x,£ )Nx= - J log0 1/n (x,E o )dx = Ai, n (E o ) 
n Jt n Jt 

> A!(£o)>k + A2(£ ) 

> k - <5 + A 2 ,„(E ) 



= k - <5 + - ( log 02, n (x, E Q ) dx 
n Jt 



Define S„ c T as the set of x e T for which 

P 



n- 1 log||APAf(E )l|-^A !> (Eo) 



> 6o (31) 



/=i 



with p = 1 or p = 2. By Lemma |3j |S„| < e _c6 o" for n large enough. On & n 
we have 

1 

-logai /fI (x,£ ) > Ai,„(Eo) - 6o > Ai(E ) -6 
n 

1 

- loga 2/ n(x,£ ) < A 2/ n(Eo) + 2<5o < A 2 (£o) + 3<5o 
n 

whence 

a^^x, E )a M (x, Eo)" 1 > e" (K - 46o) M2<k<d (32) 
Let N > n be a multiple of n (for simplicity of notation) and write 



4 N) (E ) = n ^;t^ £ ^ 



=Jc-l 
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We wish to apply the AP to this product with Aj := A^.^Eq). In order 

to do this, we take N — le 5l "j where < 61 «c 5^ is another constant. The 
union of N/n - k shifted copies of S n has measure 

< — e - c5 l n < e -( c ^l-5i)n < e -cb\nj2 

n 

Denoting this set again by S n , we see that conditions I© hold with fi deter- 
mined by J32l) , i.e., 

where C is a large constant, say C = k/(26i). By the LDT applied to Af n \Eo) 
and A<?\E ), we may also assume that dlOl l holds on & n for all < <k. 
Hence, ((TT]) holds, viz., 



k-2 k-2 



log||Af(Eo)ll + £log||^. n JEo)ll-£log||A^(E )ll < C— (33) 

;=1 ;=0 "Vf 1 

provided x e & n . The right-hand side of ((33|) is of the form N~ c for a large 
constant C, as is the measure of &„. Averaging <(33j) over S c n therefore yields 

\NA liN (E ) + {k- 2)nA lrn (E ) - 2(k - l)nA 1/2 „(E )l < N" c 

which further implies, for all n > hq, and with N = [e lJI J where 61 depends 
on k, 

|A 1/N (E ) + A ljM (E ) - 2A 1 , 2 „(E )I < = C(k)^^ (34) 



We now argue that ((34J) remains valid if we replace Eo with a nearby E. To 
see this, we first run the argument leading up to (|34l for all n € [fiQ,e n °] =: A/o 
and note that d32l in that case remains valid near Eo; in fact on some ball 
B :- B(Eq, e) c S where e depends on k and n$. For example, we may take 
£ = exp(-Ce"°) with some large C. Furthermore, the condition (O and ( fTO)) 
also remain valid on that ball. Therefore, j34l does hold for E e B with 
n e No. In particular, with 

N e {Le 6l "J I n € N ] =: M, 



we have 

logN 



A 1/N (E) > 2A 1/2 „(E) - A!,„(E) - C(k)- 



N 



logN 

>A t (E )-5 -C(K)-^- VEeB 



(35) 



The idea is now to combine this estimate with one for A 2 ,n(E). This will 
allow us to show that the gap between the {Aj in {E)\ d . =l does not shrink by 

much when we pass from scale n e No to scale N e N\ with E 6 B. If 
d — 2 we are of course done, see ((28)) . If d > 2, then we need to re-run 
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the AP-argument for A 2 A% (Eo). This is legitimate, since the two largest 
exponents of the A 2 A co-cycle are 

A 1 (E ) + A 2 (E ) > k + Ai(Eo) + A 3 (E ). 

The same AP argument as before therefore yields 

|A 2 , N (E ) + A 2/ „(E ) - 2A 2/2 „(Eo)l < = C(k)-^- (36) 

and thus, with N e N\ related to n e No as before, 

logN 



A 2/N (E) < 2A 2/2n (E) - A 2 , n (E) + C(k)- 



N 



logN 

< A 2 (E ) + 2<5 + C(k)— |— V E e B 
Subtracting this from d35t yields 

Ai / n( e )-A2 / n(E)>A 1 (Eo)-A 2 (E )-36o-C(k)^|^ VEeB (37) 
By the same token, we have this property for all gaps, i.e., for any 1 < j < d, 

A jiN (E) - A; +lfN (E) > A ; (E ) - A ;+1 (E ) - 36 - C(k)^^ > k V E e B 

provided §o is sufficiently small, and hq large. 

We claim that we may now iterate this argument without shrinking B. 
First, we define 

{Le M J \neN e ] =: N M 
for each t > 1 and note that by construction 

Uf> Wf = ["0/ °°) 

In fact, N( and JV^+i overlap. Ton & No associate N e N\ as above, and 
then set N = |_e 5lN J € N%. We may now repeat the argument leading to (|34)) , 
(|36l l, respectively, for the N, N scales. Indeed, these estimates now become 

l/WE) + A W (E) - 2A 1/2N (E)| < C(k)!^ 

} N _ (38) 

|A 2 jj(E) + A 2 , N (E) - 2A 2/2N (E)| < C(k)^|^ 

N 

for all E € B. Next, we replace Ayv/(Eo) in ((34|) (with E 6 B instead of Eo) 
with Ai /2 m(E). Subtracting the resulting estimate from ((34)) implies that 

|A 1/W (E) - A 1/2N (E)| < C(k)^ V E 6 B (39) 
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for all N € N\, and similarly for the other exponents. In combination 
with (j38|) and (|37)l this allows us to conclude that for all E € B, 

/loeN logN\ 

A^(E) - A 2 ^(E) > A 1 , N (E) - A 2/N (E) - C(k)(-|- + -§H 

N _ (40) 
JogN loeNx 
> A x E ) - A 2 (E ) - 36 - C(k)(-|— + -4r-) 

It is essential here that we do not lose more factors of b$, but only subtract 
terms such as on the right-hand side which are summable over a sequence 
of scale Ni e Nj. In view of fliPl we can indeed repeat the arguments again 

for the next scale N' :— [e 5lN ] e N3 and so on. This establishes our claim 
concerning infinite iterations while keeping the ball B fixed. 

In particular, (|39)> will hold for all large N. Summing these estimates over 
a sequence 2 k N implies that 

|A ;;N (E) - A ; -(E)| < C(k)^ V E e B (41) 

and all large N and all 1 < j < d. In particular, we have the uniform gaps 
Ay(E) - Ay+i(E) > K V E 6 B (42) 



Note that we have shown that the validity of (|42)) at one point Eq implies its 
validity for all E near Eq. 

To establish the Holder regularity of the exponents, we first replace Ai ]v(Eo) 
with Ai (E) in d34l which yields 

|Ai(E) + A 1; „(E) - 2Ai /2)3 (E)| < C- = Qk)^- V E e B (43) 

and all large N (the relation between n and N is as before). The idea is now 
to pick two E, E' £ B and then to compare the resulting estimates (|43)l with 
an N that is adjusted to t>(E,E'), the distance between E,E' 6 S. The be 
specific, we claim the crude bound 

|Ai,„(E) - Ax, n (E')| < e c "" b(E, E'f\ (44) 

for some ySo £ (0,1]. This follows from the assumed Holder regularity of 
E 1 — > A(x, E) with a a constant Co that depends on A. Indeed, expanding the 
co-cycle product and writing the difference in a telescoping fashion leads 
to the estimate 

|l4 n) (E) - 4 n) (E')ll < e C " HE, E'f° (45) 

Furthermore, 

\\A { :\E)\\>c(d)\detAf\E)\^>ce- Cn 
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whence 

1 r i \\A ( i ) (E)-A ( i ) (E , )\U 
\A h n(E)-M,n(E')\<- log 1+ " x \ n) x m )dx 
"Jt v \\A ( "\e)\\ ; 

< sup l|A ^ )(£) ~^ i)(£,)l1 < Ce Cn b(E, E't 



, 6 t \\Af{E)\\ 

for all n > 1 as claimed. 

We infer from d43b, (HU) that 



|Ai(£) - Ai(E')\ < e c ° n t>(E, E'f* + C(k)^ 



N 



< b(E, E'f° N^ 1 + C{k) 1 ^j^- VE,E'£B 



(46) 



Optimizing over N yields the desired result: 

|Ai(E) - Ai(E')l < C( K )b(E, E'f V E,E' e B (47) 

where y = y(K,fio). The same argument applies to the other exponents by 
repeating these considerations for the powers A^A^\e). □ 

It is natural to ask whether A\(E) remains Holder if Ai(E) - Ai{E) > 0. In 
other words, if there is a gap between the two largest Lyapunov exponents, 
does it follow that the top exponent is Holder? The proof we just gave 
does not show this, since it relies on statements about A2 /f! (E) over different 
scales n, uniformly in E. Such control can only be obtained via the AP, 

at least within the confines of our methods. However, the AP for A 2 A^ 
requires a gap between Ai and A3, and so on. So it is really important for 
our argument that all exponents are distinct, even to show that the top one 
is Holder regular. 

We remark that the proof we just gave is not optimal in some ways. The 
main improvements one can make relate to various upper bounds we used, 
such as ((44)) where one can insert the Lyapunov exponent Ai(E) into the 
exponent, viz. 

|A 1/n (E) - A!, n (E')| < e 1K ^ n b(E, E'f\ (48) 

for large n; here we are using that Ai(E) > which follows from J5). To 
obtain ((48)) , one relies on the upper bound 

suplog||4 n) (E)|| < nA 1 (E) + (logn) c ° 
for large n, and analogously 

v 

suplogHAMj^E)!! < nY\j(E) + (logn) c ° 
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provided that Tfj =1 A/(E) > 0. Here Co is a constant that depends on a>, and A. 
For a proof, see |GolSch2, Proposition 4.3]. Results of this type played an 
important role in the finer results required for the Cantor structure of the 
spectrum, see I GolSch2l IGolSch3| . 

Moreover, they were used by Bourgain in ||Bou2l Chapter 8] to show that 
the Holder exponent does not depend on the Lyapunov exponent (as long 
as it is positive) for Schrodinger co-cycles. In addition to these finer upper 
bounds one also needs to improve on (|2T]l to accomplish this, namely by 
removing the dependence on the 5 2 on the right-hand side. The appearance 
of this 5 2 , too, is closely related to a crude upper bound, see d26l ). Indeed, it 
is natural that we may place the Lyapunov exponent on the right-hand side 
of d26l ), which then shows that for 5 A\ we may put 5 on the right-hand 
side of ((2T|) rather than <5 2 . This is exactly what is done in [Bou2. Chapter 8]. 

However, we do not seem to gain anything from these improvements 
which is why we have not implemented them in the general GL(d, R) case. 
As already mentioned in Section [2] for d > 2 there is no relation between 
the size of the gaps and the Lyapunov exponents themselves. Thus, unlike 
the d = 2 case the methods of this paper cannot possibly lead to gap- 
independent Holder classes. 

6. Rates of convergence 

As a byproduct of the argument presented in the previous section we 
obtained the convergence rates (|4l]> . viz. 

\Aj,„(E) - A y (£)| < C(K)^p V E e £, 1 < j < d 

for all n > 2. Here £ is as in Theorem [T] (cover the compact space £ by 
finitely many open balls B as in the previous section). The purpose of this 
section is to remove the logn factor from this rate. This is the analogue 
of |GolSchlt Theorem 5.1] where the same result was obtained for SL(2, R) 
Schrodinger co-cycles. 

Proposition 2. Under the assumptions of Theorem\l\one has 

\A, n (E) - A/(E)| < - VEefi, l<j<d (49) 

" ' n 

for all n>\. The constant C depends on k,A,Q,oj. 
Proof. We write, 

Af ) (E) = A? + L(E)oA? ) (E) 
The idea is now to apply the AP to the three matrix products A^ n \E), 
A^l no) {E), and A^(E) with factors of the form A^l m w (E). Here n is large and 

Ci log n < £j < 2C\ log n, nij = ^ l\ 
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for some large C\. This can be justified as in Section [5] for all x e S c n 
where \S„\ < n~ w , say (taking C\ large). Subtracting the resulting represen- 
tation (Ell for log||A!g n(U (E)|| and log||A< n) (E)|| from that for log ||Af°(E)|| 
yields the estimate 



j log ||4 2n) (E)ll - log \\A^ n jm - log ||A< n) (E)|| 

l( 2f ) I InrrHAW fVW I _1_ 1^ 1 1 A W 



- lo S \KHn-^ + l0 § l^™(£)H + log K+(n V^ 11 I < M_1 (50) 

for all x 6 & n and some £ logn. The power n _1 can be improved on 
the right-hand side but this is of no consequence. Integrating < f50T > over T 
and noting that the integral over S n makes a negligible contribution yields 
(dropping E for simplicity) 

\2n(A li2 „ - Ai, n ) ~ 2£(A 1/2( - \ u )\ < Cn 1 
Setting R(n) := 1n\A\ i2n - A\ A \ we infer that 

C 

R(n) < R{£) + - 
n 

Iterating and summing yields R(n) < C for all ft > 1 as claimed. Applying 
the exact same reasoning to the exterior powers A? A yields the analogous 
estimate for the other exponents. □ 

We now turn to the question of whether j49l ) is an optimal estimate. As the 
example of constant A shows, this of course need not be the case. However, 
based on the observation from Section [5] that the quantities A^(E) - Ay(E) 
and Ay/(E) - Ay^E) differ only by an amount that is exponentially small in t, 
we now establish the following dichotomy: either d49|) is optimal, or the 
convergence is exponentially fast. 

Corollary 3. Under the assumptions of Theorem [I] there exists a constant c\ - 
C\{k,A,Ej,oS) > such that 

\Aj(E) - 2A m {E) + A y/ (E)| < exp(-c^) (51) 

for all sufficiently large I and all 1 < j < d. Moreover, there is l§ - £§{c\) so that if 

\A jA (E)-A j (E)\>4exp(-c 1 £ l ) 
for some t\ > £q and some j, then 

I^W £ ) - A/(E)| > ^^(E) - A ; (E)| 

for all k>0. In other words, on balls B on which (O holds, we have the following 
property: for each E e B, either Ay ;K (E) — > Ay(E) exponentially fast, or 

\A hn {E)-A ] {E)\>^- (52) 

for infinitely many n. 
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Proof. The relation ((51]) is a restatement of (|43]> . 

Now assume that A^(£) - Ay(E) > 4exp(-ci^i) where ; is fixed. In view 
of (EE), 

|A^(E) - A ; (E)| > \\A jA {E) - A ; (E)| - i exp(-crfi). 
Continuing inductively one obtains that 

\A m (E) - A ; (E)| > ^|A y/l (E) - A ; (E)| 

- 1 exp(- Cl A)[l + 2 exp(-(2 - \)c x € x ) + ... (53) 

+ 2 k ~ l exp(-(2 ,c - 1 - l)ciA)]. 
Now choose 4) so large that 

CO 

£2'exp(-(2'-l) Cl 4)<2. 

;=0 

By d53l and our assumption, 

|A //2 ^(E) - L(E)| > ^IA^(E) - A ; -(E)| - | expC-c^) 

> ^Tl A i/ 1 ( £ )- A i< E )l 
as claimed. □ 



7. Other types of shifts 

As already noted in Section |2l the proof is modular and rests on two 
main ingredients, the avalanche principle on the one hand, and the large 
deviation theorems on the other hand. In fact, the proof in Section|5]applies 
to any GL(d,K) cocycle over an "abstract" base dynamics as in Section [l] 
6 : X — > X provided one has an LDT 

V 

m({x e X | In" 1 log ||AP4' !) (E)|| - £ A;, M (E)| > S}) < exp(-c(<5)n) (54) 

i=i 

for all 5 > and all sufficiently large n > tiq(5) with some constant c(<5) > 0. 
In fact, a weaker statement such as 

V 

m({x e X | In" 1 log ||AP4' !) (£)II - £ A,>(E)| > §)) < exp(-n T ) (55) 

;'=1 

for all 6 > 0, n > fto(<5) an d some t = t(6) > would still lead to a result. 
But instead of Holder regularity one obtains a modulus of the type <[7|). For 
example, instead of j5Tb we only obtain 

|A,(E) - 2A^(E) + A^(E)| < exp(-n, < a < 1 
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under the assumption ((55)) . In the iteration which underlies the proof of 
Theorem [T] given in Section [5] we can therefore only pass from scale n to 
N — exp(n^) for some < jS < 1. The details are routine, and are left to the 
reader. 

Examples of dynamics that give rise to LDTs of the form ([55]) are shifts 
on higher tori T v as well as those on T for which ((U) is relaxed such as to 

\\nco\\ > c(co)n~ a Vn>l 

where a > 1 is arbitrary but fixed. See for example [GolSchl , Proposition 
10.2], |Bou2l Proposition 7.18] for precise statements along these lines. It 
is not known if these weaker results can be improved or not; in other 
words, for higher-dimensional tori we currently do not have an estimate 
such as (|54]> . 

8. Products of random matrices 

For the sake of completeness, and less to say anything new, we now 
sketch how the machinery of this note relates to the classical theory of 
products of random matrices. A standard reference for everything that we 
will need is the book by Bougerol, Lacroix IBouLacll ; see also Le Page [LeP], 
and Ledrappier's lecture notes IILedL 

Let jU be a probability measure on GL(d, R) and consider i.i.d. variables 
with common distribution /j such that E[log + ||Yi||] < °o. We assume 
that f.i is strongly irreducible, i.e., there is no finite union of proper subspaces 
of R d which is invariant under every matrix in Tu, the smallest closed 
semigroup containing the support of /,(. Furthermore, we assume that p. 
is contracting, i.e., there exists a sequence {Mj}j in such that ||My|| _1 My 
converges to a rank-1 matrix. 

Fiirstenberg's theorem [Fur], says that for d = 2 the co-cycle generated 
by the random sequence { Yj\ c SL(2, R) has Lyapunov exponents A\ > > 
A2 = -Ai. Moreover, there exists a unique /j-invariant probability measure 
on PR 1 , denoted by v, such that 

Ai = / f log||M-*||/i(dM)v(d*) 

JpR 1 JsL(2,R) 

where M • x is the direction of Mx where x belongs to the equivalence 
class x. If these conditions are valid for all exterior powers APY\, a theorem 
of Guivarc'h and Raugi, see [BouLac , page 78], extends Fiirstenberg's result 
to d > 2 ensuring that all exponents are distinct. 

By a theorem of Le Page [LePJ, we further know that one has exponential 
convergence to the invariant measure v, see ||BouLact page 106] for the 
precise meaning. Finally, assuming that \i has exponential moments, one 
has a large deviation estimate of the following form: for every 5 > 

limsup n~ l logP(| log \\Y„Y n -i ■ ■ ■ Y 2 Y 1 \\ -n\ 1 \> nb) < -c(6) < 

n— >oo 
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see [LeP, Theoreme 7], or |BouLad page 131]. We therefore have all in- 
gredients in order to apply the exact same argument as above. Starting 
from the case where the random co-cycle does not depend on a parame- 
ter, we conclude for example that Corollary [3] holds (with E fixed). Thus, 
one either has exponential convergence of the Lyapunov exponents, or ((52)) 
holds. However, by the aforementioned convergence result of Le Page it 
follows that ((52]) does not occur in the random case and one always has 
exponential convergence. In particular, this logic shows that the statement 
of Corollary [3]zs sharp. 

For random co-cycles depending "nicely" on a parameter, we may again 
conclude that the exponents are Holder in the parameter. However, we see 
no need to make this precise; indeed, if the dynamics is random or strongly 
mixing then we might expect much better regularity of the exponents. At 
least for Schrodinger co-cycles this is indeed the case, see Campanino, 
Klein DCamKleU and Simon, Taylor [SimT ayl- The approach we followed 
here is clearly not able to capture results of that strength; conversely, the 
techniques in these papers cannot handle deterministic dynamics such as 
shifts. 
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